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Extended Thermodynamics is a very important theory: for example, it predicts hyperbolicity, 
finite speeds of propagation waves as well as continuous dependence on initial data. Therefore, it 
constitutes a significative improvement of ordinary thermodynamics. Here its methods are applied 
to the case of an arbitrary, but fixed, number of moments. The kinetic approach has already been 
^ | developed in literature; then, the macroscopic approach is here considered and the constitutive 
functions appearing in the balance equations are determined up to whatever order with respect 



X 



to thermodynamical equilibrium. The results of the kinetic approach are a particular case of the 
present ones. 

1 Introduction 

Extended Thermodynamics is a well established and appreciated physical theory (see refs. [I], [2] 
regarding the first pioneering paper on this subject and an exhaustive description of the results which 
has been subsequently found). More recent results regarding the kinetic approach are described in 
refs [3]- [6] and many interesting properties are there obtained and exposed . The macroscopic 
approach has also been investigated but the exact solution of the condition which are present in 
the theory with many moments is still lacking. This gap is here filled and the general solution is 
obtained. 

The balance equations of extended thermodynamics with an arbitrary number of moments are 

dtF h ... in + d k F h _„ ink = Si^.An for n=0,...,N, (1) 
where we call F the tensor i^...^ when n=0. 

Here the various tensors are symmetric and F il ^ iN i : and S^...^ are supposed to be functions of the 
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previous one, in order to obtain a closed system. In particular F, Fi, Fu, Fm denote the densities 
of mass, momentum, energy, and energy flux respectively. In this way eqs. ([1]) for n = 0, 1, and the 
trace of eqs. (H|) for n=2 are the conservation laws of mass, momentum and energy; obviously to 
this end it is necessary to assume that S=0, Si = and Su = 0. 

Eq. ([1]) can be rewritten in a more compact form using a 4-dimensional notation in a space that we 
suppose to be Euclidean (nothing will change if the space is pseudo-Euclidean with -+++ signature, 
so we have chosen the simpler case). 

In particular, let us define the symmetric tensors M ai " Qjv+1 and S ai '" aN as follows: 
1. the Greek indexes go from to 3, 
2 M h...i n o...o = p h _ dn for n= o,...,N+l 

3. S*i-*»o...o = S h ... in for n=0,...,N. 
In that way the balance equations (pQ) can be simply written as 

d a M aai - aN = S ai - aN , (2) 
where d a for a = means the partial derivative with respect to time . 

The entropy principle for this equations, by using Liu's theorem [7] ensures the existence of the 
parameters L ai ... aN , called Lagrange Multipliers, such that 

dH a = L ai ... aN dM a i- a » a , L ai ... aN S a i- a » > (3) 

where H° is the entropy density and H l its flux. A brilliant Ruggeri's idea is to define 

H' a = -H a + L ai ... aN M ai " - aNa (4) 

and to take the Lagrange Multipliers as independent variables. In this way eq. ([3])i becomes 
dH' a = M ai - aNa dL ai ... aN , from which 

Q H 'a N+1 

Fit ^ ' 

"JjQl ...Ojv 

In this way the tensors appearing in the balance equations ([2]) are found as functions of the param- 
eters L ai ... aN , called also mean field, as soon as H a is known. Obviously L ai ___ aN is symmetric. By 
substituting ([5]) into eq. ([2]) this takes the symmetric form 



dL/3 v -p N dL ai . 



or r<cn---ajv 

-Oa N+ i-L'l3i-p N — & , 



■a N 



so that hyperbolicity is ensured provided that H a is a convex function of the mean field. By 
eliminating these parameters from eqs. © we obtain Fi 1 ^^ N+l again, as function of F, Fi, 
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Fi 1 ,,,i N . If we want a model in which some among eqs. (pQ) is present only by means of one of its 

traces, it can be obtained from the present model with the method of the subsystems [2]. 

Note that eq. © for aia2---«7V+i = ii---i n i n +iO...O an d for a±a2---aN+i = i\...i n Q..M n+ i gives 

respectively 

p. .. - dH '° F .. (6) 

as in the 3-dimensional notation. 

So, to impose eq. © we have to find the more general expression of H a such that M ai ° 2 ' ■■ aN + 1 is 
symmetric. We will refer to this as "the symmetry condition". We will impose also the principle 
of galilean invariance; this has been exploited in [5].[8]-[I0] for a generic system of balance laws. In 
section 2 we will apply these results to our system, taking care of converting them in the present 
4-dimensional notation, so obtaining further conditions. 

In section 3 these, together with the symmetry condition, will be investigated and their solution 
will be found up to whatever order with respect to thermodynamical equilibrium, except for two 
numerable families of constants arising from integration. In section 4 it will be shown that the results 
of the kinetic approach are a particular case of the present one so that, as usual, the macroscopic 
approach is more general than the kinetic one. If we rewrite the present paper with N-l instead 
of N, we find the model with less moments with the method used in this work, which we call the 
"direct method". But in ref. [2] it has been shown that a solution (for the model with N-l instead 
of N) can be obtained also with the "method of subsystems" ; this consists in taking the costitutive 
functions of the model with N as maximum order of moments, and in calculating them in Aj 1 ...j Ar= o, 
i.e., for zero value of the 3-dimensional Lagrange multiplier with the greatest order. In section 5 
we will see that the solution obtained with the method of subsystems is a particular one of that 
obtained with the direct method; more explicitly, it can be obtained from the latter by considering 
equal to zero one of the above mentioned families of constants. At last, conclusions will be drown. 



2 The Galilean relativity principle. 

To impose this principle, it is firstly necessary to know how our variables transform under a change 
of galileanly equivalent frames £ and This problem has been studied by Ruggeri in [8] and we 
have only to write its results in our 4-dimensional form. This is easily achieved in the kinetic model 
because the kinetic counterpart of M ai "' aN+1 is 

M ai-a N+1 = j ^ . . . c <*N+i d£ ( 7 ) 
with c° = 1, dc = dc^d(?d(? and f is the distribution function. Consequently, in £' we have 



m 



i = J f c ' ai ■■■c aN+1 dc' 
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and, if v % is the constant velocity of each point of £' with respect to E, we have c a = c a + v a , with 
v° = 0. It follows that 



N+l / 

M oi...«iv + i = I 

8=0 V 



iV + l 



or 

nr i 1 / 

N + l 



N+l / 

M ai - Qjv+1 = V] I 

i=0 V 



with ta = (1,0,0,0) for our previous notation. We obtain the transformation of M ai '" aN (which 
was the initial independent variable) multiplying eq. ([8]) by t aN+1 so finding 

M ai...a*o = X^l"-^ (v)m^- pN ° (9) 

with 



N 

yati...a N 
A 3 1 .../3 iv 

8=0 



where we have taken into account of v° = 0, of the identity I J 1 4 = I I and that the 

V * / w 

term with i=N+l gives a null contribution. Comparison between ([9]) and (llOj) with (jSj) shows that 
^^i"'^n c0iua - be obtained from simply replacing N+l with N. From eq. (|10p it follows 

also 

Similarly, if" transforms according to the rule 

H a = h°v a + h a , (12) 

of which H° = h° is a component. 

Eqs. ([8]) and (|12j) have been obtained with the kinetic model only for the sake of simplicity; it 
is obvious that they hold also in the macroscopic case. The transformation rule of the Lagrange 
multipliers can be obtained now from ([3])i with a = 0, i.e. 

dh° = dH° = L av .. aN dM a ^° = L ai ... aN X^ym^-^° 

where (|12p and ([9]) have been used. In other words we have 

dh° = lf h ...p Jf dmfr~f > »° (13) 

with 



L a\...a N — ^ v ai...a N ±J Pi---PN 
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i.e. 



N 



N 



.a N 



J2[ . v-^ ?jft - , ' ftL «.+i...^)A...ft 



(14) 



i=0 



A consequence of this result can be obtained from ([4]) with a = and written in the frame i.e 
h'° = -h° + l ai ... aN m ai - aN °; it follows dh'° = m ai - aN °dl ai ... aN from which 

dh'° 



m 



dl ai ... aN 

as in S. Moreover, from (gj), (Jl2j), ©, ((TTJ), dHJ) and again gj) and fll2j) it follows 



H 



h°i) a — h a -L T V" 1 ' 

.h°v a -h a + lp 1 ...p N m^-^\ a + l0 1 ...p N m 0l "- 0rfa 



(15) 



i.e. 



(16) 



which is similar to (|12p . 

We are now ready to consider the Galilean relativity principle. It imposes that the following diagram 
is commutative 




/ — Y ai -- aN L 



■v-Pi-Pn+i S 1 —d N+1 ( ya\...aN r ^ 
V a t S h'\---) + h' a {---) 



m 8i-6 N+1 (vai-a N t \ 

III ^.-^...^jy. U ai ... aN J 



In other words, we must have 

H' a (L ai ... aN ) = «%/i'*(x£v:;?/^ 

Eq. (JTT])2, by using eqs. (fTTj) and ([15]) becomes 



x Pi...p N dh n<a 



L 5i...5j\ 



v a + X, 



L <5i ...<5jv 

/3i— I^N^Si—Sna 



Si...8f 



m 



(17) 
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Now the derivative of (I17j) i with respect to L^'"^ N is 



It follows that eq. (HZ]) 2, holds iff 

dh 



"71— 7JV 

i.e. 



Ql 7i--'7iV 71— 7JV 



ci 7l ... 7jv 



(18) 



which is the counterpart of eq. ([5]) in the frame 

There remain to impose eq. (I17p i . Now it becomes an identity when calculated in v = (see eqs. 
(fT6j) and ([10]) to this regard) so that it holds iff its derivative with respect to Vj is satisfied, i.e., 

dh'° dU 



for a = 0, 



= / l °^ + -f T!" ,7iV for a = 0,1, 2, 3. (19) 



*'7i---7iV 

dlyi— 7jv 

^ <9/ 71 ... 7iV 

The second of this has been obtained by taking into account also eq. (|19pi ; on the other hand, this 
is included in (fT9"]) 9 with a = 0. Eq. (TIT)]) ?, by using eq. (fTljb now becomes 

<9/i' Q N f N \ 

h ' 0s ? + Wi E • * ' *( ai -i ai ^-^" 1 £ ai+1 ... aN )ft...ft- 1 j = 0. 

We remove the symmetrization with respect to a\ ■ ■ ■ aj\r which is not necessary because of the 



dh 

then reintroduce the symmetrization with respect to a\...aj^^i, obtaining so 



contraction with which is symmetric; for the same reason we can exchange and ocn and 



dh' a A / n 



h'°Sf + t aN - Yl ■ •■<(«i-V/-« ft - I Wi)A..W = °- 



ui<<xi...a. N i=l \ I 

We replace i with i+1 and we have 



dh' a ^ I N 



h '° 5 j+wi wE . n (^ + i)-^ 1 -^^ i -^ i ^ +1 ... ajv „ 1 )/5 1 .../3u = >° r 

ol ai...a N i=Q \ l + 1 y 



^U 1 „ / N- 1 
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But, by using eq. (fl4|) we have 




because tj = 0. This allows to rewrite eq. ([2"Uj) as 

dh' a 

= h%Sf + iV^ t aN l ai ... aN _ d . (22) 

Until now we have obtained that the entropy principle jointly with the galilean relativity principle 
amounts to say that 

1. eqs. (J5]) are invariant under changes of galileanly equivalent observers (see eq. (|18p ). 

2. the further condition (|22|) must hold. 

For the sake of completeness, we note that eq. (|17P i might be satisfied also with H a and h a , i.e. 

H a {La 1 --a N ) = v a tsh s (X"*,.."£ L ai -a N ) + h a {X^...^ L ai ... aN ) . 

But this is a consequence of (|17p as it can be seen running over backwards the above passages 
which allowed to obtain eq. (|16p from eq. (|12p . Moreover, in [9] and [TO] it has be proved that the 
conditions here obtained are the same of the following approach: 

1. consider eqs. ©, (fT2|) and (fT6|) but with vi = y-, instead of an arbitrary constant Vi] in 
this way m ai '" aN+1 , /i a and h a become the non-convective parts of M ai '" aN+1 , H a and H a , 
respectively, 

2. impose the conditions (fTHj) and (p2|) but considering ^ 7l ... 7JV independent variables, 

3. consider eqs. (fTHj) with a = and m 4 °'"° as definition of Z 7l ... 7JV = Z 7l ... 7JV (m ai '" Qjv °), and 
substitute this in the expressions of rn n "" lN + 1 , h a and h' a so obtaining the closure in terms 
of the non-convective quantities m a i - a N0_ 

In any case, we have to impose (|18p and (|22p ; in other words we have to find the quadrivector 
h' aN + 1 such that the right hand side of eq. (|18p is symmetric and for which eq. (|22p holds; after 
that eq. (|18p gives m^ 1 '"^ N ^ N+1 . In this way we will find the required closure satisfying the entropy 
principle and that of galilean relativity. This will be done in the next section. 
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3 Exploitation of the conditions (1181) and (1221). 



We want now to impose eqs. (|18p and (|22p up to whatever order with respect to thermodynamical 
equilibrium. This is defined as the state where 

1 

lpi-p N = Xt Pi '" t PN + s^uhpiP^fo ' ' ' 1 Pn) ( 23 ) 
holds, with hp~f = 5^ — tpt^ = diag(0, 1, 1, 1), 

A = • ••i^-Av A« = ^ /^ 2 i* • ••#%... /9w . (24) 

We can consider the Taylor expansion for 

with 

la. a.. = la. a.. - Xta. ■ --ta.. - 

3 



kV 

k=0 

hi-pN = 1 Pi-Pn - Xt Pi ■ ■ ■ t PN ~ o^^(Aft^3 • • • Hn), ( 26 ) 



f) k h' a 



- l*^s:J e , (27) 

where the multi-index notation Bi = (3} ■ ■ ■ (if has been used. Thanks to eq. (|18|) we can exchange 
a with each other index taken from those included in any Bi. So it is possible to exchange every 
index with all the others, i.e., A aBl " Bk is symmetric with respect to any couple of indexes. We 



note that there are 2 compatibility conditions between eqs. (j25[) and ()2T[) : they can be obtained as 
follows: let us consider the tensor Fi P kh y — as function of Ir.X, A//, and take the derivatives with 
respect to lp v ..p N , calculating the result at equilibrium; we find 

A aB 1 ...B h p x .-p N _ ( d k+1 h a ^ \ dl 7l ... 7jv + 
\dl Bl ■ ■ ■ dl Bk dl 11 ... lN ) eq dl^...p N 

^ ( d k+1 h' Q \ dX + ( d k+1 h' a \ dX u 



dl Bl ■ ■ ■ dl Bk dX J ^ dlp v ..j3 N ydl Bl • • • dl Bk dX a J ^ dlp v ..p N 
If we multiply this by tp 1 ■ ■ ■ tp N and by h^^t^ ■ ■ ■ tp N we find, respectively 

A aB 1 -B k (3 1 -p Nt ^ ... t p N = d_A aBl - B >' 
AaB 1 -B k f3 1 ---/3 N i, f ...+„ — Q 9 AaBi—B k 

h n P\Pl l Pz l PN— 6 d\ii ' 

where we have taken into account that from eqs. (|24p and (|26p it follows 
dX ffii... t P N dX u = ( N ) h{PiP* t fo ■ ■ ■ t M 



(28) 



dlpi-PN 91 Pv-Pn \2 / 

dL h-lN _ Q9i J3 N ) f Pi f P Nf . ../ _l N \ u(PiP2 f P 3 ... f P N ) h/ f ... f . 
ui Pi-Pn j \ 2 I 
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from which 



9A 3A 

™ • ■ • */3jv = 1 hpifotfo ■ ■ ■ t j3 N - 

d\i dXu 

■Kj */3i ■ • • tp N = — hpifotpz ■ ■ ■ tp N = 3 

u <"yi— 7jv . . _ n t/t 7l ... 7jv 



n, "*/3i ' • • *ftv - o7 h PiP2 t P 3 ' - - *ftv = °" 

° 1 Pi~Pn ° L Pl"PN 

It will be useful in the sequel to note a consequence of the condition (|28p . By using also eq. ([25 
we have 



E 1 4"£i---B fc _ 1 7i---7 Jv r 



_fPl ... fPN f _ . . . t _ _ 

3 V 2 



a (Pl ...„Pn) _ f Pl ...fPNf ...f _ I ( N | hiPlPlffh ...+M h , f ...f . 

y-yi 1 U L 1\ b lN n „ "(7172 ^73 l 7iv) 



fc=0 ' v 7 

+ £ h. {w, AaB "" B> ) r * • " k****** ■ ■ ■ ^ (3 

00 1 

E 1 AaBi-Bk^Pi-PN] . 7 
,^0^1)! 1 



where conditions (|28j) have been used in the last passage. So we have proved that derivation of eq. 
([25]) with respect to lp 1 ...p N is equivalent to its derivation with respect to lp 1 ...p N , but considering 
independent the components of this tensor, except for the symmetry. Proceeding with the subse- 
quent derivatives and calculating the result at equilibrium, we find eq. (|27p . In other words we can 
forget eq. (f2T|) but we have to retain eqs. ([28]) . We have then to transform eqs. (fTBj) . ([22]) and ([28]) 
in conditions for the tensor A aB i— B k^ the above mentioned symmetry of this tensor ensures that 
eq. (fT8|) is satisfied. Before imposing eqs. ([22]) and ([28]) . we note that the most general expression 
for a symmetric tensor depending on the scalars A, An and on t a is 

r ivfc+1 1 

L 2 J f Nk + l\ 

A a ia 2-a Nk+1 = ^ Ufc,2 S (A,A n )/i (Qia2 . . . h a2s - ia2s t a2s+1 ...^ fe +i) (29) 



s=0 



2s 



where the binomial factor has been introduced for later convenience. Thanks to this, eqs. ([28 
become 



9k+l,2s — §\9k,2s 



^9k+l,2s+2 - 2T+3 6 W^9k,2s- 



for 8 = 0, • • • , [^±±] (30) 



3 EXPLOITATION OF THE CONDITIONS ([IB) AND ([220 . 



10 



There remains to consider eq. (122|1 : thanks to eq. (125j) . (|23h and (129p . its value at equilibrium is 

2 

= 5o,o + g A «5i,2 



which, thanks to eq. (|30|) 9. becomes 

_ 8 

= 50,0 + o *u~^r- 9o,o- 
3 <9A/j 

Its solution is 



SO,o = A a f G ,o(A) ) (31) 
with Go,o(A) an arbitrary single variable function. 

But eq. (|22p is equivalent to its value at equilibrium, and to its r th derivatives with respect to l Bi 
calculated at equilibrium, for all values of r. The r th derivatives of eq. (|22p with respect to l Bi is 

cfh'^i 8 r+1 h' a 
rF/V a 81 

-|- iv/t QA , „ n n _ , yoz,) 

Ola 1 -a N Cll'(B 1 ' ' ' 8l.Br-! 8L Br) 

where the indicated symmetrization is treated as the multi-index B{ was a single index. The eq. 
(|32p can be easily proved with the iterative procedure. Now we have to calculate this expression at 
equilibrium. Let us evaluate each single term of this relation. 

• Thanks to eqs. (|25p and (|23p . we have for the first term 

ga ( d r /t% \ _ ,a AflBl -.B r . 

*> {di Bl ---di Br - 6 > A V 

\ / eq 

• The second term at equilibrium, thanks to eq. (I23p . is 



at / / I at AaBi—B r ai—<XN± _\ _u j. j 

dl Bl ■■■dl Br dl ai ... aN aN tai " a »-u J " iVA W3^K JV -^(« 1 «a 2 '"W-i) 

/ eg 

The symmetrization in the right hand side can be omitted because the term is contracted with a 
symmetric tensor. Now we use eq. ()29p . We see that the terms containing the factor t ai gives zero 
contribute, so that the above expression can be written as 

[iV(r + l) + l | 



V g r +i2s ( iV ( r + 1 ) + 1 ] 2s — h a ^ a2 ---h a2s - ia2 n a2s + 1 ---t a N(r+Dt a h aN 

^ y + ' I 2s / N(r + 1) + 1 w 



s=l 

1 

3 



' A//2/ij ftl ■ ■ ■ t ctN _ 1 



where the indexes in i?i • • • S r and a at are included into the oij; after the contraction with t a2 ■ ■ ■ t aN 
this expression becomes 
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where the indexes 7' represent -Bi ■ ■ ■ B r . 

• Let us evaluate now the contribute of the last term in eq. (|32p . i.e. 



iVrt Qjv 



•o-N-13 



dl ai ... aN dl Bl ■ ■ • 0J Br _i / ^ 



Nrt aN A aBl "- Br - iai "- aN g^ • ~ gSHrZlh^ = Nrt aN A aaNB ^- B r-^--- p *-ihf N) 



where we have esplicitated B r = j3\ ■ ■ ■ (3 r N . We can now prove that 

d r h' a \ dl 



Nrt aN 



dl ai ... aN dl( Bl ■ ■ • 3/_e r _ 1 I dl Br ) 



eq 

is symmetric with respect to two generic indexes (5f and /?*, with s < t = 1, • • ■ , r. In fact it can be 
written as 



■OiNj 



,_ l ° N dl ai ... aN dl Bl ■ ■ ■ dls^dls^ ■ ■ ■ dl Br dB k 
fc=l,--- ,r 

+ ^ t ^ J 4a«ivBi---B s -iB a +i---Bt-i^---/3^Bt + i---Br(/3f---/3| >r _ 1 / l / 3 l >r ) + 
+i Y f ^^aa w B 1 .--B s _ 1 ^--- / 3| r B s+r --S t _ 1 B t+r --B r (/3*--- / 9^_ 1 ^)_ 

The first of these terms is clearly symmetric with respect to 0? and /?*, while the sum of the last 
two is 



+terms like /„ x ,l ,u ' x r - V'7/' ; 
that is obviously symmetric with respect to /3f and /3^. 

Consequently our tensor is symmetric with respect to every couple of indexes taken between 
B\ ■ ■ ■ B r , so that it can be expressed as 

[JVrl 

/ N r\ 

Nrt aif A aarT ^"' p ^"^'"^hf I,) = 2s g r ,2sh a{12 ■ ■ ■ h^-^f 3 ^ 1 . . . t y Nr . 



s=Q 



2s 



( Nrl 

h lNr+l) + y ( Nr _ 2s) M gr 2s fc h {l2lB . . . /l 72 s 72 s+ l t 72 3+ 2 . . . f (33) 

Here we have calculated firstly t aN A aaN ^ 1 "'^ N "'^ 1 "'^ N - 1 by using eq. (J22D and then distinguishing 
the terms in which a is index of an h" from those in which it is an index of a t; finally we have 
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multiplied the result times /iJ JVr+1 and symmetrized with respect to 72 • • • 7^7-4-1 ■ 



Until now we have finished to evaluate the three terms of eq. f|32|) calculated at equilibrium; so it 
becomes 



= ^2 9r+l,2s 

8=1 



iVr + l\ 2 
.2s-l/3 



A«/lJ 72 • • • ^2,-172^72^ + 1 . . . t 7Nr+l t a) + 



+ X/ #r+l,2s 

8=1 



jVV I 1 \ 9 

-X u hY 2 • • • ^72,-172^72,+! . . . tfttrr+lf) + 

2s-lh J 



Nr 
2s . 



+ £2 

8=0 

m 

+ ^ (Nr - 2s 



9r 2 S h a{12 ■ ■ ■ ^2,-172, t 72, + l . . . + 7iVr/j7iVr + l) 



+ 



8 = 



I Nr ] 



Nr 
. 2s . 



9r,2. 



,£ a /j(7273 . . . /j72s72s + l £72s+2 . . . £JNr f{l Nr + 1 ' 



8 = 

+ £ 

8=0 



( Nr + 1) ( ^ ) 9r,2sh {al2 ■ ■ ■ h^-^ 2 H^ +1 ■ ■ ■ t ^h] Nr+l) + 



^ + 1 ) l^l9r+l,2s+2h^ 2 ■ ■ ■ ^ S+ l72 s +2 t 72 s +3 . . . ftor+lf*) 



(34) 



where in the second term we have changed the summation index s according to s=S+l. 

Note that this equation is automatically symmetric. In [TT] was proved that = is an identity 

for the case of 13 moments; here we find that this property is valid also for an arbitrary number of 

moments. 

So we have proved that eq. ([32]) amounts to 



= (Nr + 1) 



Nr 
2s , 



9r,2s + 



Nr + l\ 2 
,25 + 1/3 



hi9r+\,2s+2 , i-e., 



2 1 

9r,2s + 3 A " 2 S + l 9r+1 ' 2s+2 = fOT S = ' ' ' ' ' 



Nr 



(35) 



Consequently, all our conditions are equivalent to the scalar eqs. ([30]) . ([31]) and (135]) which are 
constraints on the scalars g r ^ s of the expansion (|29p . It remains to exploit them. For s = 0, • • • , [^f] 
we can substitute g r +± 2 S +2 from eq. ([30]) into eq. ([35]) 9 which now becomes 



2 8 

-Ai/77T— 5r,28 + 5r,28 = 



2s + 3 "9A H - 



whose solution is 



2s+3 



9r,2s = \i 2 G ri 2s(A) for s = 0, 



(36) 
(37) 
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In this way eq. (j30"]) 9 is exausted, except for s = Afr 2 +1 but only for the case with Nr odd. 

If Nr is even eq. (|37p holds for all g r 2 S , while if A^r is odd the validity of eq. (|37p is not still proved 

for <7 r> jVr+l- But for AV odd we can use eqs. (j30|) with k = r, s = Nr + l , i.e., 



9 

Q\9r,Nr+l — 9r+l,Nr+l 
8„ _ iVr+4 1 „ 



In the right hand sides we can use eq. (|37p because £^±1 < 
except for the trivial cases N=l,2. In this way the system ([38 



N(r+1) 
2 



and ^±3 < 



AT(r+l) 
2 



becomes 



_ JVr + 4 

^\9r,Nr+l = \i 2 G r +i i iVr+l(A) 

_ JVr + 6 

9 _ 7Vr+4 1 \ 2 n l\\ 



(38) 
hold, 

(39) 



The integrability conditions for this system gives 

G>+l,i\rr+3 = ^-(Nr + 2)G r+l , Nr+1 . (40) 
After that the system (|39p can be integrated and gives 

JVr + 4 

9r,Nr+l = 2 G r ,Nr+lW + C r ,Nr+l, (41) 

with 

2 1 

G r ,Nr+l = -nT~ T ~ n^Wl ,Nr+3, (42) 

3 JVr + 2 



while c rj Ar r+ i is an arbitrary constant arising from integration. So eq. (j37j) is a valid solution also 
in the case Nr odd and s = Nr + 1 : except to add the arbitrary constant c^jvy+i- 
Now we can see that this constant doesn't occur in eq. (|3(jp i (because the right hand side is 
differentiated, while in the left hand side and in the case N(k+1) odd, we have 2s < 2 f^ 1 ] from 
which 2s < N(k + 1) + 1). Nor it occurs in eqs. flS), (EH), &, (02} and ([35]) (the proof for this 
last equation amounts to verify that [%f] < for Nr odd and [%f] + 1 < [ JV(r + 1)+1 l for 

N(r+1) odd; obviously, in both of them we have N odd. If r is odd too, we have to verify only the 
first one, i.e. Nr ^~ 1 < 2 +1 , which is an identity; if r is even, we have to verify only the second one, 
i.e. ^ + K Nir+ 2 1)+1 which is true, at least for N > 1). 

On the other hand, the contribute of this constant to the tensor A aia2 "' aNk + 1 is /i( Ql ° 2 • ■ ■ h aN k aN k+^ ■ 
Ck JVfc+l) as it can be seen from eq. (|29j) . The contribute of all these constants to h a follows from 
eq. ([25]) and reads 

oo 1 

E 7^c 2r+lM2r+1)+1 h^l ■ • • ■ ■ ■ • • • • ■ • • (43) 

where we have put k = 2r + 1. 

It is easy to verify that this additional term satisfies identically the symmetry conditions for eq. 
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(I18j) and (|22p (in fact i Qjv is contracted with an h aN , for this additional term). In other words, we 
can assume eq. ([37]) for all g r ^ s (also for s = [ 2 ] )> except that, in the case with N odd, we have 
to add to h' a the additional term (|43p . 

Let's then substitute from eq. ([371) into eq. (|30]) i and (|35|) : so they become 



G 



Gk+i,28 = G'k,2s for s = 0, • • • , 
„2s + 1 



r+l,2s+2 



G r o s for s = 0, 



(44) 



(45) 



But this last equation holds also for s = 0, • • ■ , [ j ; this is obvious when Nr is even, while it is 
just eq. (|42p when Nr is odd (remember that we have eq. (|42p only for the case with Nr odd). 
After that, we see that eq. (f3Tj) is contained in ([37]) for r=s=0, while eq. ([4*0]) . by using eq. (|4*2|) . 
becomes G' N Nr+1 = GV+i,jVr+l which is just eq. (|H|) with k=r and s = [ — r 2 +1 ] (remember that eq. 
POj) holds only for Nr odd). 

There remain eqs. ([44p and (|45p . To this end, let us define H r<s from 



G 



r,2s 



(2f)! 
2 s s! 



In this way eqs. (|44p and (|4"5]) become 



r+l,s+l — H rat 



H' r ,s = -^-H r+1)S for s = 0, 



Nr + 1 



(46) 



(47) 



Eq. (|47pi suggests to define H r s also for s > . In fact, let h be a number such that 

; we can define H rjS = H r+ h,s+h- In this way eq. 



s + h < 



N(r+h)+l 



2s-Nr+l 



N-2 



(for example, h 
{\47h holds for all r and s. Regarding eq. (JITJta we have 

-3 



H'r,s — Hr+h,s+h ~ - ^ _ -^r+/i+l,s+/i — -Z~Hr+l,s', 



in other word, also (|47p ? holds for all r and s. 
After that, 



if r > s we have 



if r < s we have 



r—s,0 



Hr.K — Ht 



d r ~ s H { 



0,0 



r,s — "0,s— r- 



(48) 



(49) 



In this way H r<s is known except for Hq p . 

On the other hand, it is easy to see that ([15]) and ([4^1) satisfy eq. (f4Tl) i . Regarding (f4Tl) 9. we see 
that 
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• if r > s =4> r + 1 > s, we have to use eq. (138]) for both sides of eq. (1371) 9 and it becomes an 
identity, 

• if r = s — 1, we have to use eq. (|3U|) for the left hand side of eq. (|37j) 9 and eq. (USD for the 
right hand side. The result is Hq 1 = ^Ho t o, 

• if r < s — 1, we have to use eq. (|39|) for both sides of eq. (|37]> 9 which becomes H' Q s _ r = 

In conclusion, i7o,o is arbitrary and i?o,p is defined by 

_3 

Hq,p = -^Hq^-i, (50) 

except for a constant arising from integration, after that, eq. (|48p and (|49p give all the other 
functions H TjS . 



4 The kinetic approach 

Let us now search a solution, for conditions (|18p and (|22p . of the form 

/;' / F{lp 1 ...p N c^---c^c a dl 



(51) 



where F is an arbitrary single variable function; it is related to the distribution function, but this 
relation doesn't affect the following considerations, so that we choose to omit it. 
The symmetry for the left hand side of eq. (|18p is certainly ensured; remembering that c = 1, eq. 
(1221) becomes 







_d_ r 

d7 3 



Fllp v ..p N cfr-..c'P»)c<* 



dc' 



(52) 



which is certainly true. The expansion of eq. (|5ip with respect to equilibrium is 

k=0 J 



c' a c B ^ 



c Bk dd_l Bl ---l Bk 



where eqs. (I24p . (I26p and the multi-index notation have been used. Then we have obtained eq. (125 
with 



A 



Qk gaBi—B k 

dXl. 



D 



aB v -B k 



F 



X + lxac' 2 



■' a c B '---c Bk dc'. 



(53) 



It is easy to verify that eqs. (I28p are satisfied with this expression. The integral in eq. (153^ 9 can be 
calculated with a well known procedure. To reach faster the result, let us consider the tensor 



F 



X H A;;C' 

3 



12 



PI Ps — 
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The above tensor depends only on scalar quantities and is symmetric, so it is equal to 

if s is odd, 

g s (X, X u )h^ 2 ■ ■ • fc7—i7.) if s i s even. 

To know <7 S (A, A^) it suffices to multiply both members by /i 7l72 • • • h- fs _ 1 - ys obtaining 



A + -X u c 



12 



Js+2 



sm 



277 



dc' = g s (X,X u )(s + 1) 



where we have changed the integration variables according to the rule 



c 1 = c sin 9 cos < 



We obtain 



with 



d G [0,+oo[, 

9s(X,Xu) : 
G S (X) 



c = c sm v sm ( 
0€ [0,tt], 



c = c cos I 



G [0, 2vr[ 



4tt 

8+1 

4tt 



X + -Xud 



12 



's+2jl 



dc' = (A„)— a" G S (X) 



s + 1 

For the sequel it will be useful to note that 



rj s+2 dr], 7] = \/%ic' . 



G' S {X) 



An 
s + 1 



4tt- 



F 



drj 



rfdr} 



■q s+1 drj- 



= —{s-l)G s ^ 

provided that Fr/ S+1 is infinitesimal for rj going to infinity. After that, we have 

£71-> = B (hfk~f>r (7,71 + t ^ + ^ . . . + tft ^) 

= £ I") B*~*h% ■ ■ ■ hJ M+1 t^ ■ ■ ■ ttfrt = 



s=0 



2q + 3 



g=0 



, 2 1. 

This allows to rewrite eq. (I53p as 

rfejv+ii 

afc r>aB 1 -B k L 2 J 
q=0 



^ /l (7172 • • • . . . t >) Ajj 2 G 2q {X). 



dX k 



kN + 1 



2q + :-i 



(54) 



(55) 



^(7172 . . . ^723-172^72,+! . . . t ^t a h u 2 G { 2 k q >(X). 



This result confirms eq. (|29l) also in the kinetic case, but with gk,2s{X, Xu) = X u 2 G 2s (A), and it 



is easy to see that these functions gt 2s satisfy eq. (|30|) . as consequence of eq. (|55|h Also eqs. ([3 
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and (|41j) are confirmed, with GV,2s(A) = G^J (X), c r ,Nr+i = 0. In this way we see that the additional 
term (|43p is not present in the kinetic approach. Moreover, the matrix H TtS defined in eq. (|46p 
becomes, in this approach 



rsl G^(X) 
J2^ &2s (A) ' 



and eq. (f50l) becomes a consequence of eq. (i55j) . But the constants arising from integration of eq. 
(|50p are not present in the kinetic approach, because all the functions G S (X) are defined by (|54p in 
terms of the single variable function F. 



5 On subsystems. 

We aim to obtain now the model with N-l instead of N through the method of subsystems. To 
this end we firstly need the relation between the 4-dimensional Lagrange multipliers and the 3- 
dimensional ones. The first of these are defined by eq. ([3|), from which we obtain 

dti - L ai ... aN dM N -L N dM N g ai p x ■ ■ ■ g aN /3 N - 

= ° N dM^f 1 t3N (h ai p 1 + ta-^tfa) ■ ■ ■ {h aN /3 N + t aN t(3 N ) = 



N 



El \ J a l— a r a r+l—a N , ] i /r a/3 1 —l3 r l3 r+1 --f3 N r , 
I j l N aM N n- ai f 3l ---h ar p r t ar+1 ---t aN t l 3 r+1 ---t[3 N 

r=0 \ r / 

N 



r=0 



w 

and = M^-^-^h% ■ ■ ■ h%t Pr+1 ...tp N . (56) 

Eq. (|56[) i gives the 3-dimensional Lagrange multipliers in terms of the 4-dimensional ones. We have 
introduced the index N to remember that we are considering the model with N as maximum order 
of moments. In this way it will be distinguished from that with N-l instead of N. 
The inverse of eq. ([56"j) i is 

7C*i— a N _ ; ft-fo «i . . . «ji _ ift-feaai j. +«u \ . . . a«iv I ^ _ 

N — L N y/3i 9/3 N ~ L N V'^T" 1 z /3i) y^N^ 1 1 Pn> — 

= E ( N ) lfr~^ N h^ ■ ■ ■ hft^ ■ ■ ■ t^t Pa+1 ...t PN = 

s=0 \ S / 

N 

= Yj X N V " astas+1 ■■■t aN) . (57) 



s=0 
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The model with N-l instead of N can be obtained as subsystem of the above one by taking 

\%~ a » = 0, 

A$- a « = A^lT for s = 0, • • • , N - 1. (58) 

We have now to express these relations in terms of the 4-dimensional Lagrange multipliers; to this 
end we see that 

N-l 

Q- a " = ^ A^p^ 1 • • • t a »\ (59) 

s=0 

while eq. ([56]) i . with N-l instead of N, is 



31— jr ~ I I AT— 1 ' t Q r -Jr l 'Q: T .+i b a N _ 1 - \W) 



Then, by substituting eq. ([511}) in eq. (|STJj) we find 

N-l 

jai-a N _ ST^ Acts+i +a Nh ai La a )/7i-7a7s+i-7JV-l + ^ 

'AT — Z-^ 71 ' ' ' 7s N-l b ls+i ' ' ' l 1n-i 

from which 



N-l 



l N ~ l N-l 

because 



l%- a « = Z^™- 1 ^), (61) 



] a 1 ---a N _ 1 ,a N _ pi-lN-l j.a N /lcii , ±ai± \ . . . ( f 1 °' N ~ 1 -L t a N-if \- 
l N-l L ~ l N-l b \""yi ' 1 l li) \""1n-i > b l 1n-\) ~ 

N - 1 fN-\\ 

ST~ \ I i7i""7s—7JV-l/ J (ai . . . h^sf fOs+i . . ,f f a N _x) f a N 

— \ \ N-i Si 'Ss Ts+i' nn-i 1 1 



s=0 

Now, from eq. ([23]) . we have 



3 



This and eq. (1231) yield 



tai---ajv /( Q l-" a JV-l j.ajv) 

eg. — eg. E ' 



that is, eq. (|6ip holds also when we calculate it at equilibrium. The deviation of eq. (|6ip from its 
value at equilibrium is 

Jai-a N = j(a 1 ^a N . ltaN ). ^ 

in other words, eq. (|6ip holds when we substitute the Lagrange multipliers with their deviation with 
respect to equilibrium. We can now substitute eq. (|62p into eq. (|25p ; in this way we find the coun- 
terpart of (|25|) with N-l instead of N. To this end we have to contract an index of each Bi, ■ ■ ■ ,B^ 
with at.; in other words, we have to contract the expression ([29]) with t a(N _ 1)k+2 ■ ■ ■ t aNk+1 . It is easy 



6 CONCLUSIONS. 



19 



to verify that in this way eq. (129p remains unchanged except that now N-l replaces N; obviously 
this is true also for gt,2s where s now goes from to ( Af ~i) fc + 1 _ This property is transferred to 
Gk,2s for eq. ([37|) and to H r>s for eq. (JM|). But H rjS is defined by eqs. ([4*8]) and (f4U|) in terms of 
Hq }P which are determined by eq. (|50|) , Therefore, the family of constants arising by integrating eq. 
([50]) . is inherited also by the subsystem. 

We have only to notice that from eq. (|49p it follows that Hq^ p is useful for H r ^ r+P which, for eq. (|46p 
is useful for G r ^( r +p)- It follows that Hq jP is present in the subsystem when r+p < 



(jV-l)r+l 
2 



is p < 



(jV-3)r+l 
2 



(jV-2)r+l 
2 



that 

Now, for a fixed 



while for the initial system was useful when p < 
value of p, it is always possible to find r such that both of the previous inequalities are satisfied. 
The only difference is that in the subsystem, Hq^ p occurs only in terms of higher order with respect 
to equilibrium, than in the initial system. This is true, provided that N > 3, that is if neighter the 
system, nor the subsystem are the 10 moments model. 

But what happens to the other family of constants, that is for the supplementary term (|43p ? 
If N is even, the model has not this term and, consequently, it cannot be inherited by the subsystem. 
If N is odd, this term is present; but when we substitute eq. (|6ip in (|43p we obtain zero because 
each i. is contracted with a projector h". We expected this result because, with N odd, we have 
N-l even in which case the term (|43p is not present. We may conclude that the other family of 
constants, or the supplementary term (|43f) . disappears in the subsystem. Only the other family of 
constants is inherited. 

This can be seen also from the following viewpoint: the family of constants arising from integration 
of eq. ([50]) . in the case N=3, will perpetuate also for the subsequent values of N; equivalently, we 
can say that the closure in the model with a generic N > 3 is exactly determined in terms of that 
with N=3, except for the supplementary term 



6 Conclusions. 

In this paper we have applied the method of Extended thermodynamics to the case of an arbitrary 
but fixed number of moments. This case has already been developed in the kinetic approach. Here 
we have considered the macroscopic approach and the costitutive functions have been determined 
up to whatever order with respect to thermodynamical equilibrium. In this way we have been able 
to find the exact general solution of the conditions present in extended thermodynamics for the case 
of many moments. The results founded in the kinetic case, as we expected to found, are a particular 
case of the present ones. Moreover, we have introduced an innovative 4-dimensional notation that 
simpliflies very much the form of the equations. Finally, the present results in the case N=3 are 
equivalent to those found in |12j . 
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